Implication of the PAMELA antiproton data for dark matter indirect detection at LHC 
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Since the PAMELA results on the "anomalously" high positron fraction and the lack of antiproton excess in 
our Galaxy, there has been a tremendous number of studies advocating new types of dark matter, with larger cou- 
plings to electrons than to quarks. This raises the question of the production of dark matter particles (and heavy 
associated coloured states) at LHC. Here, we explore a very simple benchmark dark matter model and show 
that, in spite of the agreement between the PAMELA antiproton measurements and the expected astrophysical 
secondary background, there is room for large couplings of a WIMP candidate to heavy quarks. Contrary to 
what could have been naively anticipated, the PAMELA p/p measurements do not challenge dark matter model 
building, as far as the quark sector is concerned. A quarkophillic species is therefore not forbidden.Owing to 
these large couplings, one would expect that a new production channel opens up at the LHC, through quark- 
quark and quark-gluon interactions. Alas, when the PDF of the quark is taken into account, prospects for a 
copious production fade away. 



I. INTRODUCTION 

The new results of experiments such as PAMELA OHl], ATIC 
H, HESS fl, FERMI H have attracted a lot of attention in 
the dark matter community during the last few months @- 
§| and somehow revived hopes that cosmic ray data could fi- 
nally shed light on the nature of dark matter. In particular, the 
apparent conflict between the seeming excess in the positron 
fraction and the lack of visible anomalies in the antiproton 
spectrum has encouraged theoreticians to propose new dark 
matter models with very unusual properties. For example, to 
fit PAMELA data, a large amount of dark matter models with 
very large values of the annihilation cross section into leptons 
and small values of the annihilation cross section into quarks 
have been proposed in the literature. All these models rely on 
"local" enhancement mechanisms so that the candidate abun- 
dance remains compatible with the measured dark matter relic 
density while, simultaneously, leading to a very high amount 
of positrons in our galaxy. The complexity of the mechanisms 
involved illustrates the difficulty to fit the data and reflects the 
ongoing theoretical effort to explain the PAMELA positron 
excess. However, this positron anomaly could be the signature 
of nearby pulsars supplementing the interstellar medium with 
a radiation of electrons and positrons. The possibility that 
spallation reactions take place during the acceleration of pri- 
mary cosmic rays has also been suggested. In this Letter, we 
concentrate on the antiproton measurements which are subject 
to neither speculation nor controversial interpretation. It has 
been shown indeed that the PAMELA p/p ratio is compatible 
with a pure astrophysical explanation. Antiprotons are gener- 
ated inside the galactic disc and result from the interaction of 
high energy cosmic ray protons and helium nuclei on the inter- 
stellar material. The lack of excess in antiprotons may suggest 
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that the WIMP couplings to quarks must be small. However, 
here, we propose to find what is the maximal coupling that 
is allowed by the PAMELA antiproton data to verify whether 
this assertion is true or not. 

For simplicity, we base our analysis on a generic dark mat- 
ter model where the WIMP % is directly coupled to a Stan- 
dard Model quark q and a heavy coloured partner F q . Inspired 
by supersy mine try, we shall focus this Letter on a Majorana 
WIMP so F q is in fact a scalar (and would be equivalent to a 
squark in supersymmetry). Our approach is generic enough to 
be applied to other types of dark matter. However a complete 
survey of the various possibilities (including scalar dark mat- 
ter species % and fermionic partner F q ) is beyond the scope 
of this work (although we anticipate that the results presented 
here are fairly general). 

As we will show, the couplings that we obtain are in fact large 
enough to also have to address the question of the implications 
of the PAMELA antiproton data for the prospects of dark mat- 
ter indirect detection at LHC. 

In "Standard" supersymmetry, the heavy states F q are pro- 
duced mostly through gluon-gluon fusion proton-proton colli- 
sions. They are expected decay into missing energy (the dark 
matter) and a jet (corresponding to the quark q) shortly af- 
ter production. Hence, this channel constitutes an important 
source for dark matter particles production at LHC. 
However, if the direct dark matter couplings to quarks are 
larger than usually expected in "standard" supersymmetry, 
new channels could open up with qq, qq and qg interactions, 
thus providing new possible signatures and increasing the dis- 
covery potential at LHC. Hence the importance of character- 
ising the direct dark matter couplings to quarks in light of 
PAMELA antiproton data. 

To constrain the dark matter characteristics (mass and pro- 
perties of the particles to which dark matter is coupled), we 
proceed as follows: we first use in Sec.|ll]a realistic cosmic 
ray propagation code to calculate the antiproton flux at the 
earth expected from conventional spallations and dark matter 
annihilations. The latter depend on the % — q — F q left and 
right couplings cf and c q R which we want to constrain from 
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the PAMELA antiproton measurements. We then consider a 
variety of annihilation channels and derive the maximal cross 
sections allowed by the PAMELA antiproton data. Results are 
expressed in units of the canonical thermal value of 3 x 10~ 26 
cm 3 s _I and referred to as the boost factor hereafter. Because 
couplings to u, d and s quarks are already severely constrained 
by direct detection, we focus on c and b quarks. With the help 
of a Monte Carlo Markov Chain, we delineate in Sec. HITIthe 
region of parameter space which saturates the upper bound on 
the boost. Finally, in Sec. [IVj we determine the number of 
events associated with qq — > FF, qq —> FF* and qg — > F% 
production at LHC and discuss the propects for detecting dark 
matter related signatures at LHC. We conclude in Sec. [V] 

H. THE PAMELA ANTIPROTON CONSTRAINTS 

The antiproton flux at the earth arises from two sources. Con- 
ventional (secondary) antiprotons are produced by the inter- 
actions of high energy cosmic ray protons and helium nuclei 
with the gas of the galactic disc. The annihilation of hypothet- 
ical dark matter species in the halo of the Milky Way generates 
additional (primary) antiprotons. To compute both si gna ls, we 
have considered the same propagation model as in lUOIl with 
the parameters best fitting the boron to carbon ratio fiUl . As 
regards the secondary component, we have used the same lo- 
cal proton and helium nuclei fluxes as in |[l2ll together with 
the radial distribution of SN remnants given by IILJl to retro- 
propagate these fluxes all over the diffusive halo. 
Concerning the primaries, we have considered a dark 
matter halo computed by lfl4ll with a local density of 
P0 = 0.3 GeVcm~ 3 . The distance from the earth to the 
galactic center is 8.5 kpc. The annihilation cross section is set 
to the conventional thermal value mentioned above. We have 
determined the boost by which that value can be increased 
without exceeding the PAMELA data. We have scanned the 
WIMP mass from 100 GeV to 1 TeV and considered dif- 
ferent annihilation channels. For each channel, the antipro- 
ton spectrum before propagation has been calculated with the 
PYTHIA 6.4 program 111511 . The results are summed up in 
Tab.U 
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TABLE I: Maximum annihilation cross section allowed by the 
PAMELA antiproton measurements for various WIMP masses and 
different annihilation channels. This cross section is expressed in 
units of the canonical value of 3 x 10~ 26 cm 3 s _1 . 



These constraints are quite severe. However, antiproton flux 
calculations suffer from a lot of uncertainties and changing 
some of the choices we made may affect these results. As 
recalled in fl6Tl the local value of the dark matter density is not 
very well constrained and in fact lies in [0.2; 0.9] GeV cm~ 3 . 
The primary antiproton flux being proportional to the square 
of the local density, one can divide all the results of Tab. Uby 
a factor ranging from ~ 0.4 to 9. The other uncertainties are 
summed up in Tab. [TT] 
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TABLE II: Maximum boost (i.e. annihilation cross-section in units 
of the canonical thermal value of 3 x 10~ 26 cm 3 s _1 ) allowed by 
the PAMELA antiproton measurements for various WIMP masses 
in the case of annihilation into bb pairs when varying the parame- 
ters. Cases 1 and 2 correspond to extreme propagation parameters in 
agreement with B/C. For case 2 the boost has been rounded to the 
closest decade. Cases 3, 4 and 5 correspond to various dark matter 
halo profiles : NFW (J3] (3), Moore (H (4) and an isothermal cored 
profile (5). Cases 6 and 7 correspond to alternative fits of the injec- 
tion proton and helium nuclei spectra respectively proposed by flOll 
and tH. 

One should notice that for most cases, the constraints come 
from the point of highest energy (61.2 GeV) published by 
PAMELA. This point is the one that suffers the biggest statis- 
tical error and its systematic error is unknown yet. Hence the 
results may change with future publication of new PAMELA 
data. Indeed the current data correspond to only 500 days of 
data collection starting the 15th of June 2006 but the satel- 
lite is still in orbit and should carry on taking data for at least 
few more months. Moreover, the uncertainty related to the 
injection spectra should diminish as soon as absolute fluxes 
for protons and antiprotons are published. Finally, new data 
are also expected from PAMELA for the boron to carbon ra- 
tio that should help us limit the uncertainties on propagation 
parameters. 



III. THE MCMC AND THE ALLOWED REGION IN 
PARAMETER SPACE 

To find the region of parameter space that reproduces the boost 
factor as determined in the previous section, we perform a 
Markov Chain Monte Carlo search. The free parameters of 
our benchmark model are a priori the dark matter mass, the 
masses of the heavy scalar partners F q , and their couplings c[ 
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and c R to Standard Model quarks q. Since b-quarks are a bit 
easier to tag (and direct detection experiments indicate that the 
couplings to u, d and s quarks must be suppressed), we shall 
only focus on predictions associated with b quarks. The only 
free parameters of interest are therefore the mass otk of the b - 
coloured state Fb, the dark matter mass m% and its couplings to 
the b quark, namely c b L and c R . Our purpose is to explore this 
4-D parameter space in order to delineate the regions where 
the boost is close to the maximal value allowed by PAMELA. 
The expression of the annihilation cross section is borrowed 
from 12011 . The range of masses and couplings which we con- 
sider lies between [100,2500] GeV and [0.01,3] respectively. 
A larger upper limit on the couplings would induce a loss of 
perturbativity. 

Translating the maximal boost B max found in the previous sec- 
tion into regions of the parameter space is tricky. For this 
purpose, we have built a Monte Carlo Markov Chain which 
explores the parameter space and affects to each of its points 
a likelihood depending on whether or not the boost is close to 
the maximal value allowed by PAMELA. This method allows 
to rapidly delineate the interesting regions. The likelihood is 
defined as a Gaussian centered on B max with width \>b = 0.01. 
Results are displayed in the correlation plots of Fig. Q] The 
first four rows correspond to the parameters m z , nif b , c b L and 
c b R of the model while the fifth one is the reconstruction of the 
boost factor. Notice that the first plot of this row corresponds 
(as it should) to the column bb of Tab. Q] 

The panel in the very first line represents the dark matter mass 
while the last panel of the second, third and fourth rows stands 
for nif b , c b L and c b R respectively. The first panel of the second 
line is a correlation plot between mp b and m x . The first (sec- 
ond) plot in the third line features the correlation between c b L 
and m-f_ ( c h L and nif b ) and so on. The dotted lines in all the 
plots feature the prior distribution while the solid line is the 
posterior distribution, i.e., the values of the parameters that are 
found by the Markov chain. The chain has explored enough 
points since the prior distributions are indeed rather flat for the 
four parameters while it matches a Gaussian distribution for 
the boost. We can also check that the contraint of dark matter 
stability (nif b > m x ) is correctly reproduced in the (m^m^) 
plot of the second row. 

Points of the parameter space are selected according to how 
the boost is close to B max . The MCMC considers a priori 
equally light and heavy WIMPs. But the heavier the WIMP, 
the larger the value of B max as is clear in Tab. U as well as in 
the first plot of the fifth row of Fig. Q] On the other hand, a 
heavy WIMP means an even heavier coloured partner Fb and 
since the annihilation cross section is a decreasing function 
of mp b , a large boost is only recovered for very large values 
of the couplings c b L and c R . Hence, we expect the MCMC 
to select large values for these couplings. This is confirmed 
by our results since both left and right handed couplings can 
reach two or three units (with 90% and 68% CL respectively). 
If we consider a larger variance \'b as in Fig. [2] the trend is 
the same. Although slightly smaller values of the boost B can 
be in principle achieved with simultaneously large values of 
the couplings and mass nif b , the latter would have to be larger 



than the 2.5 TeV limit which we set. Hence, in this case, our 
Monte Carlo chooses a large value of one of the coupling (e.g. 
c£) together with a smaller value of the other coupling (e.g. 
c R ) and a large value of nip. (still within our limits). Our main 
conclusion is that the PAMELA antiproton observations do 
not preclude large values for c b L and c R . Our naive benchmark 
model provides therefore a simple example where the WIMP 
can be quarkophilic and yet satisfy the PAMELA constraints. 
Note that the correlation plots displayed in Fig. [T] are ob- 
tained by saturating the upper limit on the boost obtained 
in section [TTJ The parameters found in these plots there- 
fore correspond to a large annihilation cross section into 
quarks, which translates into a relic density smaller than 
a dm h 2 = 0.1099 ± 0.0062 JH. Since our model simply re- 
lies on WIMP couplings to quarks (we do not invoke the Som- 
merfeld mechanism for example nor annihilations into heavy 
gauge bosons), this means that we have delineated the region 
of the parameter space where WIMPs have the "wrong" relic 
density. Discovering heavy coloured states at LHC decaying 
into jet + missing energy and with properties matching this 
region of the parameter space would therefore jeopardise con- 
ventional dark matter scenarios. 



IV. HEAVY STATES PRODUCTION AT LHC 

Since the WIMP couplings to b quarks can be very large, the 
question of F q production at LHC from quark-quark collisions 
arises. We are interested in fact in bb — s- FbFb (which may hap- 
pen if dark matter is a real scalar or, in our case, a Majorana 
particle), bb — > FbFb and qg — > F%. Interestingly enough, the 
qq — > FF and qg — > F% cross sections can reach up a few pb 
or a few hundred fb respectively for large couplings. 
Such large cross sections seem rather encouraging. However, 
to get realistic estimates, one has to take into account the par- 
ton distribution function (pdf) associated with b-quarks and 
gluons. By focusing on b, we have definitely considered a 
case where there is a very large suppression due to the pdf 
l22tl . However, the pdf for c quarks is also suppressed and 
the boost factor is not so large. Hence b quarks are more or 
less representative of what is to be at most expected at LHC 
for these types of Fb production processes, given PAMELA 
antiproton measurements. 

At large x and Q 2 , convoluting our cross section with the 
parton distribution function leads to a large suppression of 
the number of events. We predict less than ten events for 
pp — > FbFb (for a integrated luminosity of fb^ 1 ) in the re- 
gion where the bb —> FbFb process is supposed to be maximal 
(i.e. for large couplings) 112411 . The pp — > F% cross section is 
a bit less suppressed because it involves gluon. However, it is 
still very small (it can reach at most ten fb, assuming a very 
small uncertainty on the value of the boost factor). Assuming 
a greater uncertainty on the boost value does not help. The two 
cross sections become even smaller since, in this case, larger 
nif b values are preferred when the dark matter couplings to 
quarks is greater than unity. Both cases lead to a too small 
number of events to be detected via the decay (jet+missing 
energy) of the heavy coloured states Fb- Indeed, they would 



4 



correspond to pure hadronic final states which are extremely 
difficult to exploit. 



V. CONCLUSION 

In this letter, we have investigated the constraints that 
PAMELA set on the WIMP couplings to quarks. We found 
that coupling values larger than the strong coupling con- 
stant could be allowed by PAMELA antiproton measure- 
ments. This encouraged us to look at the production of heavy 
coloured final states from proton proton collisions. We fo- 
cus in particular on collisions producing Ft particles (and 
therefore b + %) in the final state, since they might be easi- 
er to tag than c-quarks. Also we ignore couplings to u,d,s 
quarks since dark matter direct detection constraints indicate 
that these must be suppressed. 

We find that, despite the large values of the WIMPs couplings 
to quarks, there is very little hope to produce a large number 
of coloured states Fy through bb, bb,bg process in pp colli- 
sions because of the b quark pdf suppression (less than a few 
100 events at LHC for most of the parameter space under con- 
sideration and a luminosity of lfb~ l ). Also, such a final state 
would be purely hadronic and would be difficult to disentan- 
gle from background. Hence, the best channel to constrain 
the WIMP (direct) couplings to quarks in the absence of other 
signatures may be, indeed, heavy coloured states production 
via gluon fusion and their decay into jets plus missing energy. 
Although there have been many studies advocating leptophilic 
dark matter particles to fit the PAMELA positron excess, our 
analysis shows that the PAMELA antiproton measurement 
also allows for "quarkphilic" dark matter in principle. This 



property may hold whether dark matter is at the origin of the 
positron excess or not. To fit both PAMELA antiproton data 
and positron excess, it therefore seems likely that only the 
WIMP couplings to gauge boson (the W in particular) may 
have to be suppressed. Not only are large WIMPs couplings 
to heavy quarks allowed by PAMELA antiproton data but they 
may also be difficult to constrain using LHC data since their 
effect on the cross section can be compensated by very large 
F q masses and is suppressed by heavy quark pdf folding. 
This actually constitutes the key point of our analysis and sug- 
gests that the dark matter production channel in this simple 
model will also rely on F q production through gluon gluon fu- 
sion. Hence, such "quarphilic" models may be difficult to dis- 
entangle from standard supersymmetry, unless one can mea- 
sure the couplings very accurately. In any case, this analysis 
shows that even though leptophilic dark matter particles may 
fit the PAMELA positron excess data, this does not imply that 
the dark matter couplings to quarks must be suppressed. On 
the other hand, the strength of the dark matter-quark interac- 
tions will remain relatively weak owing to the exchange of 
heavy messengers F q . 

Acknowledgments 

We are grateful to Peter Skands who was very kind to help 
us in dealing with Pythia code, J. Hamann and J. Lesgour- 
gues for their help with getdist routine. We would like to 
thank F. Arleo, F. Boudjema, J. Idarraga and W. Porod for 
fruitful discussions. T.D. is grateful for financial support from 
the Rhone- Alpes region (Explora'Doc) and from the Interna- 
tional Doctoral school in AstroParticle Physics (IDAPP). F.S. 
is supported by the DFG Graduiertenkolleg GRK-1 147 while 
the work of R.K.S. is supported by the German BMBF under 
contract 05HT6WWA. We would like to thank the referee for 
his/her suggestions. 



[1] O. Adriani, G. C. Barbarino, G. A. Bazilevskaya, R. Bellotti, 
M. Boezio, E. A. Bogomolov, L. Bonechi, M. Bongi, V. Bon- 
vicini, S. Bottai, et al., Nature (London) 458, 607 (2009), 
0810.4995. 

[2] O. Adriani, G. C. Barbarino, G. A. Bazilevskaya, R. Bellotti, 
M. Boezio, E. A. Bogomolov, L. Bonechi, M. Bongi, V. Bon- 
vicini, S. Bottai, et al., Physical Review Letters 102, 051101 
(2009), 0810.4994. 

[3] J. Chang, J. H. Adams, H. S. Ahn, G. L. Bashindzhagyan, 
M. Christl, O. Ganel, T. G. Guzik, J. Isbert, K. C. Kim, E. N. 
Kuznetsov, et al., Nature (London) 456, 362 (2008). 

[4] H. E. S. S. Collaboration: F. Aharonian, ArXiv e-prints (2009), 
0905.0105. 

[5] A. A. Abdo, M. Ackermann, M. Ajello, W. B. Atwood, M. Ax- 
elsson, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, M. Bat- 
telino, et al., Physical Review Letters 102, 181101 (2009), 
0905.0025. 

[6] L. Bergstrom, T. Bringmann, and J. Edsjo (2008), 0808.3725. 
[7] D. Hooper, A. Stebbins, and K. M. Zurek (2008), 0812.3202. 
[8] M. Lattanzi and J. I. Silk (2008), 0812.0360. 
[9] M. Cirelli, M. Kadastik, M. Raidal, and A. Strumia, Nuclear 
Physics B 813, 1 (2009), 0809.2409. 



[10] F. Donato, D. Maurin, P. Salati, A. Barrau, G. Boudoul, 
and R. Taillet, Astrophys. J. 563, 172 (2001), arXiv:astro- 
ph/0103150. 

[11] D. Maurin, F. Donato, R. Taillet, and P. Salati, Astrophys. J. 
555, 585 (2001). 

[12] F. Donato, D. Maurin, P. Brun, T. Delahaye, and P. Salati, Phys. 

Rev. Lett. 102, 071301 (2009), 0810.5292. 
[13] W. R. Webber, M. A. Lee, and M. Gupta, Astrophys. J. 390, 96 

(1992). 

[14] J. Diemand et al. (2008), 0805.1244. 

[15] T. Sjostrand, S. Mrenna, and P. Skands, JHEP 05, 026 (2006), 
hep-ph/0603175. 

[16] G. Jungman, M. Kamionkowski, and K. Griest, Physics Reports 

267, 195 (1996), arXiv:hep-ph/9506380. 
[17] J. F. Navarro, C. S. Frenk, and S. D. M. White, Astrophys. J. 

490, 493 (1997), astro-ph/9611107. 
[18] J. Diemand, B. Moore, and J. Stadel, Mon. Not. Roy. Astron. 

Soc. 353, 624 (2004), astro-ph/0402267. 
[19] Y. Shikaze et al., Astropart. Phys. 28, 154 (2007), astro- 

ph/0611388. 

[20] C. Boehm and P. Fayet, Nucl. Phys. B683, 219 (2004), hep- 
ph/0305261. 



5 



[21] J. Dunkley et al. (WMAP), Astrophys. J. Suppl. 180, 306 

(2009), 0803.0586. 
[22] M. R. Whalley, D. Bourilkov, and R. C. Group (2005), hep- 

ph/0508110. 

[23] A. Lewis and S. Bridle, Phys. Rev. D66, 103511 (2002), astro- 



ph/0205436. 

[24] G. Belanger et al, Comput. Phys. Commun. 177, 894-895 
(2007). 



6 




200 600 1000 
1600 r— ■ ■ — 




200 600 1000 500 1000 1500 




200 600 1000 500 1000 1500 1 2 3 




200 600 1000 500 1000 1500 1 2 3 1 2 




200 600 1000 500 1000 1500 1 2 3 1 2 3 5 10 15 20 




200 600 1000 500 1000 1500 1 2 3 1 2 3 5 10 15 20 -2 -1 1 




200 600 1000 500 1000 1500 1 2 3 1 2 3 5 10 15 20 -2 -1 1 -1 1 2 

m ? [GeV] i*i f [GbV] c j> ^ B tog 10 [a(pp -> F D FJ/fb] log |Q [o(pp -» F^fc] 



FIG. 1: The dark matter parameter space is tunned to reproduce the boost factor obtained in Sec.[II|for a bb pair. The central boost value is 
taken from Tab.[T]while its uncertainty is taken to be vb =0.1. Due to this very small boost uncertainty, one can see the very strong correlation 
between the boost value and the dark matter mass (see Tab.0. The cross sections account for the F^F^, F^F^, F^F^ and F^F/y as well as the F^x 
and Fb% final states respectively. We used the getdist routine from the COSMOMC code f23ll to make this plot. 
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FIG. 2: The dark matter parameter space is tunned to reproduce the boost factor obtained in Sec.[II|for a bb pair. The central boost value is 
taken from Tab. [J Its uncertainty is estimated from the smallest difference between the boost value displayed in the last five columns of Tab. [Hi 
and the central value of reference. Whatever the dark matter mass, the uncertainty vb thus obtained is much larger than in Fig.[T]and one loses 
the strong correlation between the value of the boost and the dark matter mass. Such "large" values of vg relaxes the constraints previously 
obtained and opens up the parameter space. As a result, large coupling values (c^ r € [1,3]) allow for larger mp b values which translate into 
smaller production cross sections at LHC. Here we use ^=10 TeV. 



